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ABSTRACT 

Temporal  and  spatial  power  spectra  from  computer  simulations  of 
the  equations  modeling  Type  II  irregularities  in  the  equatorial  elec- 
trojet are  presented.  The  power  spectra  are  computed  using  the  code 
of  Perch  and  Sudan  (1977),  averaging  over  times  as  long  as  20  seconds. 
The  dependence  of  the  average  frequencies  and  linewidths  on  wavelength 
and  angle  from  the  electrojet  drift  is  studied.  The  mean  frequencies 
obtained  agree  well  with  linear  theory  except  at  short  wavelengths 
while  the  linewidths  are  greater  than  those  predicted  by  the  linear 
damping  rates.  The  spatial  power  spectra  were  also  computed  and 
found  to  be  approximately  isotropic.  Favorable  comparisons  are  made 
with  radar  observations  and  the  theory  of  Sudan  and  Keskinen  (1977). 
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INTRODUCTION 

In  the  equatorial  Li  region,  radar  backscatter  observations  have 
shown  irregularities  in  plasma  density  whose  spectral  width  is  often 
comparable  to  tire  mean  Doppler  shift  and  whose  phase  velocities  are 
on  the  order  of  the  electron  drift  speed  [Balsley,  1969;  Balsley  and 
Farley,  1971,  1973;  Farley  and  Balsley,  1973;  Fejer  et  al.,  1975, 

1976].  Many  of  the  features  of  these  Type  II  irregularities  have 
been  explained  with  the  linear  theory  of  the  gradient  drift  insta- 
bility [Maeda  et  al.,  1963;  Knox,  1964;  Reid,  1968;  Register  and 
D’Angelo,  1970;  Whitehead,  1971].  Both  one-dimensional  [Sato  and 
Tsuda,  1967;  Sato,  1971,  1973;  Rogister,  1972]  and  two-dimensional 
[Sudan  et  al. , 1973]  nonlinear  theories  have  been  proposed  to  explain 
the  presence  of  the  Type  II  irregularities  with  wavelengtlrs  that  are 
predicted  to  be  linearly  damped.  The  two-dimensional  model  has  found 
support  both  from  experiment  [Balsley  and  Farley,  1973;  Farley  and 
Balsley,  1973;  Fejer  et  al.,  1975,  1976]  and  simulation  [McDonald 
et  al. , 1974,  1975;  Ferch  and  Sudan,  1977], 

Much  of  our  understanding  of  the  physics  of  the  Type  II  irregu- 
larities is  based  upon  VI1F  radar  observations,  particularly  the 
spectral  measurements.  However,  existing  radars  have  studied  only  a 
few  wavelengths  whereas  the  range  of  scale  sizes  of  the  Type  II  irregu- 
larities is  known  to  span  at  least  two  orders  of  magnitude.  Spectra 
from  numerical  simulations  can  be  computed  for  many  wavelengths  and, 
as  a result,  supplement  experimental  observation. 


We  have  extended  a previous  numerical  study  [Perch  and  Sudan, 
1977]  and,  using  the  same  code,  computed  the  temporal  and  spatial 
power  spectra  from  simulations  of  the  two-dimensional  model.  This 
code  is  sufficiently  fast  that  moderate  spectral  integration  times 
can  be  achieved.  The  numerically  obtained  temporal  power  spectra 
favorably  reproduces  Type  II  radar  spectra. 

In  this  paper  we  first  review  the  linear  and  nonlinear  theory 
of  the  Type  II  irregularities.  Then  we  present  the  results  of  the 
numerical  simulations  while  making  comparisons  with  experiment  and 
theory.  Finally  we  summarize  and  discuss  the  principal  results  of 
this  study. 


LINEAR  THEORY 

For  wavelengths  of  more  than  a few  meters,  the  E region  electrojet 
plasma  can  be  described  as  an  electron  and  singly  charged  ion  fluid, 

I The  fluid  equation  can  be  simplified  by  making  the  following  assump- 

tions : 

1.  The  electron  gyroradius  is  small  compared  to  wavelengths 
of  interest. 

2.  The  wave  electric  field  is  essentially  electrostatic,  since 
the  plasma  pressure  is  much  smaller  than  the  ambient  magnetic  energy 
density. 

3.  The  electron  and  ion  densities  are  equal  (quasineutrality) 
since  we  are  interested  in  wavelengths  much  larger  than  the  Debye 
length. 

4.  The  ions  and  electrons  are  isothermal. 

5.  The  waves  are  assumed  to  propagate  exactly  perpendicular  to 
the  magnetic  field,  since  these  waves  suffer  the  least  diffusive 
damping. 

6.  The  electron  inertia  is  neglected,  because  the  electron- 
neutral  collision  and  cyclotron  frequencies  are  much  larger  than  the 
wave  frequency  Doppler-shifted  to  the  electrojet  drift  frame. 


7 . Ion  inertia  is  also  neglected  for  the  same  reasons  as  in  (6) . 
With  these  restrictions  the  fluid  equations  are: 
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2 Vne 


— (E  + V x B)  + u„  — + v V = 0 
m v-v*  ~e  omt  e n e*e 


(2) 


Ze 


0 Vn. 
2 i 


^ (E  + V.  x B)  + nf  — + v-V.  = 0 

m.  w -i  •«  in-  i*«a. 


(3) 


V*J  = 0 

w 


(4) 


where  a denotes  electrons  fe)  or  ions  (i)  and  J = ne(ZV--V  ),  E = -Vd> 

W vl  *45  /** 

is  the  electric  field,  B is  the  geomagnetic  field,  u,  is  the  thermal 
velocity,  and  v is  the  collision  frequency  with  neutrals.  All  other 
symbols  have  their  conventional  meaning.  We  adopt  a Cartesian 
coordinate  system  in  which  the  magnetic  field  is  aligned  with  the 
x axis  (north),  the  y axis  is  westward,  and  z direction  vertically 
upward.  Eliminating  algebraically  the  velocities  Vc,  V. , equations 
(l)-(4)  reduce  to 


+ Uax  x V4>-Vn  = DaV2n 


gu  x x V<{>*Vn  + u (V<J>*Vn  + nV2<J>)  = D’V2n 


(5) 

(6) 


where 


% = eB/me 
= ZeB/nu 
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uf.  = eft  /m  (ft2  + v2) 
H e'  ev  e e 


Dl  = vOJ2/^2  + v2) 


In  the  lower  E region  equatorial  electrojet,  the  parameters  are 
such  that  B = (1  + tJO  Vj/ft^,  ya  = (1  + i p)  1y^,  Dfl  = 2^(1  + ij>)  1D^_, 

D'  = (1  + ilO'^/ft^,  ty  = vevi/fie«i. 

Equations  (5)  and  (6)  can  be  linearized  about  the  equilibrium 
nQ(z)  » Nq(1  + z/L),  d<Po/dz  = a(z  + L)'1  where  a = B(V2/fte~ (l+(ve/fte)2)LVd] 
and  is  the  electron  drift  velocity  (in  the  y direction)  at  z = 0. 
Fluctuations  n = n exp{i[ky.y  + kzz  - (w^  + iy^)t]}  with  k*B  = 0 and 
kyL  » 1 have  frequency  and  growth  rate  given  by  [Register  and 
D'Angelo,  1970] 


“kr-=  C1 


Yk  = fod  + «[')‘1H(Vve^1VVd/k2L(1  + ^l-k2cs/vi}  (8) 


? 2 2 
where  C = K(T  + T.)/m.  = (m  /m-)n  + n.  is  the  acoustic  velocity. 

3 6 11  016  1 

Instability  is  possible  if 


V,  k , y k2C2 


Note  that  L is  positive  for  density  increasing  with  altitude  and 
Vj  is  positive  for  electrons  drifting  toward  the  west,  the  usual 
daytime  conditions. 
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NQNLINEAR  HEORY 

The  minimum  unstable  wavelength  for  the  ionospheric  parajnetcrs 
of  the  K-region  from  Fq.  (9)  is  about  30  m.  The  radar  backscatter 
observations  indicate  very  large  nonthermal  density  fluctuations  at 
3 m.  These  fluctuations  are  stable  and  damped  according  to  the 
linear  theory  of  the  previous  section  and  therefore  should  not  be 
excited.  In  addition,  vertically  propagating  irregularities  are 
observed  that  cannot  be  accounted  by  linear  theory.  The  explanation 
lies  in  the  fact  that  the  unstable  oscillations  with  wavelengths  longer 
than  30  in  couple  nonlinearly  to  these  short  wavelengths  and  in  the 
steady  state  the  rate  of  excitation  just  balances  the  damping.  How- 
ever, one  dimensional  theories  of  nonlinear  wave  coupling  [Sato  and 
Tsuda,  1967;  Tsuda  and  Sato,  1967;  Rogister,  1972]  cannot  explain 
either  the  vertical  Type  II  radar  echoes  [Farley  and  Balsley,  1973; 
Balsley  and  Farley,  1973;  Fejer  et  al.,  1976]  or  the  density  fluctuation 
amplitudes  of  2-55  observed  in  rocket  experiments  [Prakash  et  al.  1969, 
1970,  1972]  for  weak  electrojet  conditions  when  presumably  Type  II 
irregularities  dominate. 

It  was  shown  by  Sudan  et  al.  (1969)  that  a one  dimensional 
nonlinear  wave  steepening  requires  abnormally  large  fluctuating 
amplitudes  30-505  and  therefore  the  actual  mechanism  for  such  inter- 
action must  be  two-dimensional  in  nature.  Since 
the  group  velocity  of  all  the  waves  is  the  same  viz.  V^/(l  + tji)  the 
wavepackets  can  couple  strongly  to  each  other  because  they  stay  together 
for  a long  time.  In  the  frame  moving  with  velocity  V^/ (1  + \ p)  with 
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respect  to  the  laboratory  an  unstable  wave  appears  to  be  a purely 
growing  wave  since  the  real  part  of  the  frequency  doppler  shifts  to 
zero.  A horizontally  propagating  wave  will  have  local  horizontal 
density  and  potential  gradients  dh/dy  and  d$/dy.  Where  the  phase  of 
the  wave  is  such  that  (dn^Vdy)d4>^)dy  > o,  secondary  short  wavelength 
oscillation  propagating  in  the  vertical  direction  will  be  driven 
unstable  in  exact  analogy  to  the  primary  waves  driven  unstable  by 
the  equilibrium  gradients  (dno/dz)  (d4>0/dz)  > 0.  In  this  manner  it 
is  possible  to  generate  shorter  wavelength  secondary  waves  and  it  is 
also  evident  that  the  spectrum  of  these  secondary  waves  is  likely  to 
be  isotropic  in  the  plane  transverse  to  the  static  magnetic  field. 

This  theory  can  account  for  vertically  propagating  irregularities 
and  also  predicted  the  two-dimensional  turbulent  nature  of  the  lype  II 
irregularities  experimentally  verified  by  Farley  and  Balsley,  1973; 
Balsley  and  Farley,  1975  and  by  Fejer  et  al.,  1976. 

The  frequency  spectra  of  the  3 m fluctuations  measured  by  radar 
backscattering  techniques  indicate  that  the  peak  of  the  spectra  ocarrs 
at  a frequency  which  is  given  by  Eq.  (7).  However,  the  frequency  width 
of  the  spectra  is  much  larger  than  that  given  by  the  damping  rate  (8) 
and  indeed  it  is  comparable  to  This  is  also  indicative  of  the 


strong  interaction  between  waves  leading  to  a nonlinear  damping  far 
in  excess  of  linear  damping  rate  and  may  be  explained  as  follows. 
Suppose  we  launch  a test  wave  with  wave  vector  k into  a region  of  well 
developed  turbulence  containing  waves  of  all  possible  wavelengths 
propagating  in  every  direction.  The  test  wave  will  immediately  couple 
to  two  other  waves,  k',  k"  whose  wave  numbers  satisfy  k = k*  + k", 

W A/V  /w  /yy 


I 
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because  the  phase  variation  of  the  test  wave  will  he  synchronized 
to  the  vector  sum  of  the  other  two  waves.  Because  the  waves  are 
almost  dispersionless  and  the  group  velocity  of  all  the  waves  is 
approximately  equal  such  interaction  could  occur  over  a long  time 
and  the  transfer  of  energy  from  k to  k'  and  k"  would  result  in 

/w  'W 

rapid  damping  of  the  test  wave. 

Hie  quantitative  analysis  of  the  nonlinear  effects  described 

above  has  been  undertaken  by  Sudan  and  Keskinen  (1977)  using  a theory 

of  strong  turbulence  whose  basic  premise  is  the  "direct  interaction 

approximation"  pioneered  by  Kraichnan  (1959).  The  principal  results 

of  this  investigation  are:  (1)  The  steady  state  two  dimensional 

density  fluctuation  spectrum  1^  is  isotropic  and  follows  a power  law 

k"n  where  n lies  between  2 and  4,  but  is  closer  to  3.2  if  one  accepts 

radar  data  which  shows  that  Ij,  = at  3 m.  The  mean  r.m.s.  density 

fluctuation  predicted  by  this  theory  falls  within  the  range  of 

observations.  (2)  The  nonlinear  damping  rate  or  the  frequency  width 

2 2 1/2 

of  the  spectrum  is  given  by  (yj,  + Tj.)  where, 

Ty.  = const,  x Ij/^  (10) 

where  8 = (1  + if>)  and  the  constant  is  of  order  unity.  This 

expression  is  verified  against  radar  observations  and  the  simulation 
results  presented  later. 

On  the  other  hand,  Rognlien  and  Weinstock  (1973,  1974)  have 


estimated  the  saturated  amplitude  of  large  scale  gradient-drift 
fluctuations  using  the  nonlinear  theory  of  perturbed  particle  orbits 
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pv'einstock  and  Williams,  1971;  Sleeper  and  Weinstock,  1972].  The 
perturbed  orbit  theory  is  based  on  the  supposition  that  an  excited 
spectrum  of  incoherent  waves  causes  the  stochastic  dif  p_  Ion  of 
collisionless  orbits  which  leads  to  nonlinear  wave  danping.  Since 
the  ions  are  collisions!  this  concept  cannot  apply  to  them  while  the 
magnitude  of  this  effect  on  electrons  is  insignificant  because  the 
electron  gyroradius  is  negligible  compared  to  unstable  wavelengths. 
Furthermore,  such  theories  cannot  explain  the  generation  of  short 
wavelength  or  vertically  propagating  waves  since  explicit  mode 
coupling  is  not  taken  into  account. 
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NUMF.RICAL  SIMULATIONS  AND  RESULTS 

Equations  (5)  and  (6)  were  solved  numerically.  A description 
of  the  code  and  methods  of  solution  is  given  in  Ferch  and  Sudan 
(1977).  With  the  exception  of  the  drift  velocity,  V^,  identical 
parameters,  typical  of  the  daytime  equatorial  electrojet,  were  used: 

B = 3 x 10'5  W/m2,  «.  = 96  sec'1,  fie  = 5.3  x 106  sec'1,  vi  = 2.5  x 103 

-1  , in4  -1  2 nri5  2 , 2 2 c c in9  2,2 

sec  , = 4 x 10  sec  , u^  = 10  m / sec  , ug  = 5.5  x 10  m /sec  . 

Table  1 summarizes  the  models  from  which  temporal  and  spatial 
power  spectra  were  computed.  Each  model  was  initialized  with  an 
unstable  horizontal  wave  together  with  a smaller  amplitude  vertical 
perturbation.  All  models  could  be  advanced  in  time  until  a quasi- 
steady state  had  been  achieved.  Model  1,  the  most  weakly  driven, 
required  the  most  time  (approximately  32  sec)  to  saturate.  All 
three  models  were  run  on  a CDC  7600,  with  Model  1,  for  example,  taking 
about  1.5  hours  of  computer  time. 

2 

As  t increased  the  region  of  nonzero  |Sn(k),t)|  cascaded  to 
larger  values  of  (k ko,  kz/ko)  than  those  originally  excited  (1,0) 
and  (0,1);  kQ  = 2-rr/L  is  the  fundamental  wavenumber  and  L = 128 
is  the  system  size.  The  time  evolution  of  the  maximum  total  density 
fluctuation,  amplitude  of  the  initial  128  m horizontal  wave,  and  the 
maximum  horizontal  and  vertical  velocities  for  Models  1-3  are  shown 
in  Figures  la-lc.  Density  contour  plots  for  Model  3,  the  most  strongly 
driven  (Vj  = 125  m/s),  at  t = 0,  2,  4,  12,  26,  34  see  are  displayed  in 
Figures  2a-2f.  Density  contour  development  for  Models  1 and  2 is 
similar  to  Model  3. 


r * 
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The  behavior  of  Model  3 is  illustrated  in  figure  le-lf  and 
2a- 2f.  The  maximum  total  density  fluctuation  and  maximum  hori- 
zontal and  vertical  velocities,  which  peak  in  an  earlier  time 
interval,  t = 3-5  sec,  than  the  other  less  strongly  driven  models, 
attain  values  of  approximately  14.85,  425  m/s  and  300  m/s,  respec- 
tively, whereupon  they  relax  and  describe  a quasi-steady  state. 

By  t = 1 sec,  the  maximum  horizontal  and  vertical  velocities  were 
comparable  to  the  original  drift  velocity.  At  t = 2 sec,  small 
scale  turbulence  is  seen  from  the  density  contour  plot.  In  Models 
1 and  2,  similar  behavior  is  observed  but  on  a longer  time  scale 
with  smaller  density  fluctuations  and  velocities.  In  Model  1,  the 
density  fluctuations  and  horizontal  and  vertical  velocities  peak 
near  t = 6-8  sec  reaching  values  of  11.85,  250  m/s,  and  200  m/s, 
respectively. 


r 
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SPATIAL  POWER  SPECTRA 

Let  the  spatial  power  spectra  in  the  quasi-steady  state  be 

defined  as  I Qc, t)  = (L/2tt) 2< 1 6nQc , t) /nQ J 2>  where  6n  = n-nQ(z) 

2 

and  L is  the  area  of  the  mesh.  In  the  quasi-steady  state  the 
ensemble  average  is  understood  to  be  equivalent  to  a time  average. 

In  order  to  analyze  the  degree  of  isotropy  of  the  spectrum,  we 
expand  in  Fourier  series  as  follows, 

N 

I(k,t)  = I(k,t)  + l (fn(k,t)  cos  nQ  + gn(k,t)  sin  nQ) 
n=l 

where  the  sum  over  n is  restricted  to  positive  integers  and  k = |k|. 
Box  normalization  and  periodic  boundary  conditions  have  been  employed. 
The  Fourier  coefficients  are 


I(k,t)  = 


2tt 

2?  { de 


2ir 


fn(k,t) 


■H 

■ i r 


dQ  I(k,t)  cos  nQ 


gn(k»t)  = - j dQ  I(k,t)  sin  nQ 
o 

From  symmetry  considerations  (6n(y,z,t)  real)  the  odd-numbered  coef- 
ficients vanish.  The  angular  integrals  are  computed  in  discrete 
(k  ,kz)  space  using  a block- integration  nearest  grid-point  scheme. 

The  fn(k,t)  and  gn(k,t)  were  evaluated  for  k in  the  range  1 k/kQ  ^ 14. 
This  is  equivalent  to  a wavelength  regime  A = 2Tr/k  of  9 m ^ A ^ 128  m. 
Although  in  all  three  models  the  cell  spacing  was  2 m,  spectral 
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infonnation  for  ^ 9 in  was  disregarded  due  to  insufficient  spatial 

resolution  and  aliasing  limitations.  Initially,  a 15-teim  fitting 

function  was  used,  i.e.,  an  isotropic  term  and  equal  numbers  (7)  of 

sin  n0  and  cos  n0  harmonics  with  n = 2,  4,  14.  Equal  numbers 

were  used  due  to  the  suggested  isotropy  of  the  spatial  power  spectra 

[McDonald  et  al.,  1974,  1975;  Ierkic  et  al.,  1977;  Ferch  and  Sudan, 

1977] . Figure  3 shows  sample  Fourier  coefficients  computed  from 

spectra  averaged  over  2 sec  near  the  beginning  of  the  quasi -steady 

state  of  Model  2 and  indicate  an  approximate  power  law  in  k for 

1 ^ k/kQ  14.  Similar  results  are  found  for  Models  1 and  2. 

Absolute  Fourier  coefficients  a , b can  then  be  defined  using 

n n 

the  best-fit  spectral  indices,  Pn>  q^  in  the  following  manner: 

I(k,t)  = a k n + y (a  k ^n  cos  n9  + b k Si  sin  n0) 

' o L v n n 

n 

with  the  an,  bn  found  from  least  squares  fitting.  For  time  averaged 
spectra  taken  during  the  middle  and  end  of  the  quasi-steady  state 
of  Model  2,  Table  2 gives  the  best  fit  spectral  indices  and  relative 
magnitudes  (.an/aQ , bn/aQ)  the  best  fit  coefficients  that  exceeded 
0.1  and  indicates  that  the  power  spectra  are  approximately  isotropic. 
Similar  results  were  obtained  for  the  other  two  models.  The  best 
fit  spectral  indices  were  computed  using  a least  squares  fitting 
routine  with  errors  equal  to  the  standard  deviations.  The  uncertainties 
in  the  relative  coefficients  reflect  the  statistical  error  in  the 
mean  of  the  time-averaged  spectral  estimates. 

L 
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The  computed  power  laws  of  the  spatial  power  spectra  I (k)  ^ k n, 
n = 3-4,  for  A:9-128  in  are  in  agreement  with  previous  rocket  data 
(n  = 3 + 1;  A:l-15  m,  Prakash  et  al.,  1969,  1970),  recent  numerical 
studies  (n  = 3.5;  X : 9- 75  m,  McDonald  et  al.,  1974,  1975;  n = 3 + 1; 

A:9-128  m,  Perch  and  Sudan,  1977)  and  dimensional  arguments  (n  3, 

Ott  and  Farley,  1974).  In  addition  the  isotropic  nature  of  the 
computed  spatial  power  spectra  for  A: 9-128  m compares  favorably  with 
recent  radar  measurements  of  Type  II  irregularities  at  A = 3 m 
[Ierkic  et  al.,  1977]. 

It  is  of  interest  to  study  in  more  detail  the  isotropic  part  of 
the  k-spectra.  Figures  4a-4d  illustrate  the  time  evolution  of  the 
angle -averaged  power  spectra,  viz.,  I(k)  = (2tt) ^ dQ  l(k,0)  at 
t = 0,  2,  13,  25  sec  for  Model  2 (V^  = 100  m/sec).  Initially,  only 
one  wavenumber  k = 1 is  excited  (k  is  scaled  in  terms  of  the  fundamental 
kQ  = 2 /128  in  ■*■) . The  initial  power  residing  in  k = 2 to  14  is  not 
physical  but  numerical  in  origin  and  results  from  discrete  Fourier 
transformation.  In  the  linear  stages  of  the  simulation,  the  power 
in  the  initially  excited  unstable  mode  increases.  As  the  turbulence 
develops,  the  power  in  the  linearly  damped  portion  of  the  spectrum 
(k  > kc  = 3)  also  increases.  Finally,  nonlinear  interactions  in  the 
fully  developed  turbulence  has  reduced  the  power  in  the  initial  wave- 
number  and  brought  the  spectrum  roughly  in  accordance  with  a power  law. 

In  Fig.  5,  for  the  quasi-steady  states  of  Models  1-3,  I(k)  is 
plotted  as  a function  of  drift  velocity  for  several  wavelengths. 

These  curves  indicate  that  the  spectra  follow  an  approximate  power 

law  in  V^,  i.e.,  I(k)  « V^.  Unfortunately,  due  to  numerical  instabilities 
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and computational  cost,  a larger  range  of  drift  velocities  could  not 
be  studied.  From  radar  backscatter  observations  at  3 m,  a quadratic 
dependence  of  scattering  cross  section  on  drift  velocity  has  also 
been  inferred  [Balsley,  ’969;  Farley  and  Balsley , 1973], 
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TBPORAL  POWKR  SPHCTRA 


The  temporal  power  spectra  I(k,uf)  « <|Sn(k,ui)  | > were  computed 
(see  Appendix  A)  using  the  direct  method  [Blackman  and  Tukey,  1959]. 
Power  spectral  estimates  from  different  time  blocks  were  averaged 
to  reduce  the  variance.  The  sanpling  frequency  was  chosen  to  be 
always  greater  than  twice  the  frequencies  in  which  there  was  measurable 
power  in  order  to  minimize  aliasing.  In  many  cases  the  spectra  were 
further  smoothed  over  adjacent  channels  using  a Harm  window.  Final 
spectra  at  a particle  (ky,kz),  i.e.,  at  a particular  wavelength  and 
angle  to  the  electrojet  drift,  consisted  of  an  average  over  20.5  sec 
in  the  quasi- steady  states  of  the  three  models.  Several  sample 


spectra  were  computed  in  order  to  estimate  a sampling  frequency  which 

-1 


was  taken  to  be  50  sec  . Both  128  and  256  point  transforms  were 
used.  By  choosing  ry  = ky/kQ,  ry  = kz/kQ  to  lie  between  0 ^ n^  ^ 14 


and  -14  ^ nz  4:  14,  spectra  are  available  from  considerations  of 


symmetry  that  cover  the  range  0°  < 0 < 360°  and  9 m < X < 128  m where 


-1, 


X = 27T/k  and  0 tan  (k  /k  ).  Figures  6a-6c  give  sample  computed 

y z 


spectra  together  with  radar  obtained  spectra  [Balsley  and  Farley, 

1971]  taken  in  the  daytime  equatorial  electrojet  for  equivalent 
conditions.  For  X = 9 m many  of  the  computed  spectra  were  of  poor 
quality  when  integration  times  were  less  than  2 sec.  All  spectra 
vvere  normalized  to  their  maximum  values.  The  principal  characteristics 
of  the  spectra  are  their  mean  frequencies  and  average  widths.  We  have 
analyzed  for  the  three  models  the  frequencies  and  widths  of  the  computed 
spectra  throughout  the  range  of  X and  0. 
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M1AN  SPECTRUM  FREQUENCIES 

The  spectrum  frequency,  io/2tt,  is  taken  to  be  the  average  of  the 
upper  and  lower  half-power  points.  Extraneous  peaks  at  zero  fre- 
quency were  disregarded.  In  Figures  7a- 7c  are  displayed  for  Made  Is 
1-3  (Vj  = 75,  100,  125  m/s)  the  computed  spectrum  frequencies  as 
a function  of  k at  the  following  fixed  angles  to  the  electrojet 
drift  (also  the  principal  radar  angles):  90°  (vertical),  k 45°, 
i 30°.  Also  included  at  each  angle  are  the  frequencies  of  modes 
that  propagate  within  7.5°  of  these  chosen  angles.  Positive  angles 
refer  to  upward  moving  waves , negative  angles  to  downward  moving 
waves.  Because  some  spectra  were  noisy  all  points  could  not  be 
paired  at  the  opposite  angle.  The  linearly  damped  (k  > kc)  and 
growing  (k  < k ) modes  are  indicated  along  with  the  frequencies 
derived  from  the  linear  dispersion  relation.  The  maximum  error 
in  each  spectrum  frequency  is  estimated  to  be  0.2  Hz.  At  the  angles 
shown,  the  average  spectrum  frequencies  computed  from  the  quasi- 
steady turbulent  state  at  the  three  drift  velocities  over  the  range 
of  wavelengths  9-128  m are  not,  within  error  bounds,  appreciably 
different  from  the  linear  results  except  at  large  k (short  wave- 
lengths) . At  these  angles,  aside  from  the  vertical  (90°),  the  com- 
puted spectrum  frequencies  are  proportional  to  wavenuirber,  to  a k. 
Figure  8 illustrates  for  Model  2 (Vj  = 100  m/s)  the  good  agreement 
between  the  computed  phase  velocity,  to/k,  in  the  steady  turbulent 
state  and  the  predicted  dependence  from  linear  theory,  viz., 
to/k  = Vd(l  + <10  * cos9.  Similar  fits  are  found  for  the  other  two 


models. 


-1S- 


These  results  imply  that  even  in  the  quasi-steady  turbulent 
state  with  turbulent  velocity  fluctuations  6V  on  the  order  of  the 
mean  drift  (see  Fig.  la-lc),  the  average  frequencies  can  be 
found  with  the  linear  dispersion  relation  using  the  drift  velocity 


AVERAGE  SPECTRUM  WIDTHS 

The  average  spectrum  width,  Am/27r,  is  defined  as  the  difference 
between  the  upper  and  lower  half-power  points.  In  Figures  9a-9c 
are  displayed,  for  the  three  drift  velocities,  the  mean  spectral 
width  as  a function  of  k at  the  same  angles.  Linearly  growing  and 
damped  wavelengths  are  also  indicated  and  the  linear  damping  rates 
shewn.  The  maximum  error  in  each  computed  spectrum  width  is 
estimated  to  be  0.4  Hz.  At  the  angles  shown,  there  is  a definite 
trend  for  the  average  widths  to  increase  with  k,  for  all  three 
drift  velocities,  a feature  also  noted  in  radar  backscatter  obser- 
vations at  3 and  9 m [Balsley  and  Farley,  1971].  In  addition,  the 
widths  at  all  angles  and  for  all  three  drift  velocities  exceed  the 
linear  damping  rates.  Mareover,  at  fixed  k,  the  width  of  the  spectra 
of  the  linearly  clamped  wavelengths  increase  with  increasing  drift 
velocity,  a feature  not  predicted  by  linear  theory.  There  seems  to 
be  no  appreciable  difference  between  the  mean  widths  computed  for  the 
upward  moving  waves  as  compared  to  those  of  the  downward  moving 
waves  except  perhaps  at  large  k.  The  exact  dependence  of  Am  on  k, 
throughout  the  range  of  wavelengths  studied,  is  not  well-defined, 
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although,  for  example,  the  mean  widths  of  the  vertical  inodes  (90°) 
for  Model  2 fit  well  a power  law  in  k,  Aw  'v  kn  with  n = 0.8.  Figure 
10  shows  Aw/w  at  0 = 45°  for  Models  1-3.  Over  the  range  of  k, 

Aw/w  exceeds  0.5  indicating  strongly  turbulent  fluctuations.  This 
is  in  rough  agreement  with  radar  observations  [Balsley  and  Farley, 

1971],  Similar  plots  are  found  for  the  other  two  models. 

To  examine  the  dependence  of  angle- averaged  spectral  width  on 

wavenumber  k for  wavelengths  not  previously  studied  with  radar,  we 

have  plotted  in  Figure  11  the  numerically  obtained  angle-averaged 

frequency  spread  vs.  k for  the  three  drift  velocities  = 75,  100, 

125  m/sec.  At  all  three  drift  velocities,  the  angle- averaged  widths 

increase  approximately  monotonically  with  k.  We  have  fitted  these 

2 2 1/2 

data,  for  each  drift  velocity,  with  a function  of  the  form  (Yj.  + r^) 
where  Yj.  is  the  linear  growth  rate  and  is  the  predicted  linewidth 
from  Eq.  (10).  This  function  gives  a reasonable  fit  after  adjusting 
rk  by  a constant  of  order  unity.  In  this  connection  it  would  be  useful 
to  study  the  spectral  characteristics  of  wavelengths  greater  than  9 m 
by  radar. 

The  numerically  obtained  angle- averaged  spectral  width  as  a 
function  of  drift  velocity  for  several  wavelengths  is  displayed 
in  Figure  12.  At  theseewavelengths  the  data  shews  that  the  angle- 
averaged  frequency  spread  is  roughly  proportional  to  with  m = 1.5- 2. 5. 
From  radar  observations  [Balsley,  1969;  Farley  and  Falsley,  1973]  it 
has  been  inferred  that  I(k)  « giving  for  the  angle- averaged  spectral 
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broadening  Tj.  « Vj.  This  agi'ees  reasonably  well  with  the  curves 
shown  in  Figure  12.  'Hie  experimentally  observed  lincwidth  Acj  from 
9 m TYpe  II  radar  backscatter  observations  [Balsley  and  Farley, 

1971]  is  shown  in  Figure  13  as  a function  of  drift  velocity  Vj. 

The  dependence  of  Aco  with  appears  linear  as  opposed  to  the 
quadratic  results  indicated  by  theory  and  the  present  numerical 
simulations.  However,  it  must  be  noted  that  these  results  are 
obtained  from  one  zenith  angle  and  not  angle- averaged.  Furthermore, 
the  effects  on  spectral  width  of  the  changing  scattering  volume  is 
not  well-defined  throughout  the  range  of  drift  velocities  shown. 
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TEST  WAVES 

'Hie  magnitude  of  the  power  spectra  of  the  short  wavelength 
irregularities  (A  < 10  m)  is  much  smaller  than  the  fundamental 
(A  = 128  m)  by  as  much  as  four  or  five  orders  of  magnitude.  In 
order  to  study  more  closely  the  turbulent  damping  at  the  shorter 
wavelengths,  test  waves  of  the  form  A sin(ky,y  + kzz  - m^t)  were 
launched  in  the  quasi-steady  states.  The  decay  rates,  t^",  of  the 
test  waves  were  found  to  be  of  the  order  of  the  spectral  width  Am 
of  the  power  spectra  computed  at  the  wavelength  of  the  test  waves. 

The  wavelength  A and  the  propagation  angle  6 of  the  test  waves 
were  taken  to  be  9,  10  m and  ± 45°,  respectively,  with  the  amplitude 
A fixed  at  0.04.  This  amplitude  corresponds  to  the  approximate 
quasi-steady  state  amplitude  of  the  turbulent  density  fluctuations. 
Figure  14  gives  an  example  of  the  decay  of  a sample  10  m test  wave 
launched  in  the  quasi- steady  state  of  Model  2.  Table  3 compares  the 
nonlinear  damping  as  computed  from  the  spectral  width,  test  wave 
decay  rates,  and  the  theoretical  prediction  for  wavelengths  studied 
in  Models  1-3.  The  agreement  between  spectral  width  and  test  wave 
decay  rates  is  reasonable.  Better  agreement  cannot  be  expected  since 
the  spectra  were  computed  by  averaging  many  samples  whereas  the  test 
wave  decay  rates  are  found  using  only  one  sample.  However,  the  trend 
of  increasing  spectral  width  with  decreasing  wavelength  is  also  noted 
in  the  decay  rates  of  the  large  amplitude  short  wavelength  test  waves. 
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CONCLUSION 

We  have  computed  the  temporal  and  spatial  power  spectra  from 
numerical  simulation  of  Type  II  irregularities  in  the  daytime 
equatorial  electrojet.  The  overall  level  of  turbulent  density 
fluctuation  appears  to  increase  with  the  electrojet  drift  velocity, 

Vj,  in  the  range  75-125  m/s.  Appreciable  excitation  of  vertical 
and  oblique  waves  is  seen  at  drift  velocities  of  75  and  125  m/s. 

This  supplements  similar  results  of  Perch  and  Sudan  (1977)  at 
= 100  m/s  and  provides  further  support  for  the  two-dimensional 
model  of  Sudan  et  al.  (1175). 

At  the  three  drift  velocities  studied  (75,  100,  125  m/s)  and  in 
the  wavelength  range  of  approximately  9-128  m the  computed  time- 
averaged  spatial  power  spectra  I (k)  are,  on  the  average,  found  to 
be  isotropic  and  follow  a power  law  in  k(I(k)  « k n,  n = 3.5  ^ 0.5). 
This  is  in  agreement  with  experimental  results  [Prakash  et  al.,  1969, 
1971,  1972;  Ierkic  et  al.,  1977  for  X = 9 m],  numerical  simulations 
[McDonald  et  al.,  1974,  1975;  Ferch  and  Sudan,  1977]  and  dimensional 
arguments  [Ott  and  Farley,  1974].  In  addition,  some  evidence  is 
found  which  suggests  that  Ij,  “ V^,  in  agreement  with  Type  II  radar 
observations  [Balsley,  1969;  Farley  and  Balsley,  1973], 

The  numerically  computed  temporal  power  spectra  compare  favorably 
with  several  9 m Type  II  radar  spectra  taken  in  the  daytime  equatorial 
electrojet  [Balsley  and  Farley,  1971].  The  numerical  spectra  are 
computed  with  known  electro jet  parameters  with  the  radar  spectra 


being  integrated  over  a small  region  of  the  electrojet  due  to  limited 
resolution  (1-3  km).  In  addition,  only  a limited  number  of  inter- 
acting modes  can  be  kept  in  the  simulation  due  to  finite  grid  size 
and  computer  storage  limitations  and  spectral  integration  times  must 
be  minimized.  These  factors  make  direct  comparison  between  numerical 
and  radar  spectra  difficult  to  assess. 

In  the  turbulent  quasi-steady  states  of  the  three  models  studied, 
the  mean  spectrum  frequencies  of  the  numerically  obtained  spectra 
agree  well,  within  error  bounds,  with  linear  theory,  except  at  short 
wavelengths.  The  average  spectrum  widths,  on  the  other  hand,  in  the 
wavelength  range  9-128  m are  greater  than  the  linear  damping  rates 
and  increase  with  decreasing  wavelength.  Also,  at  a fixed  wavelength, 
the  average  widths  of  the  linearly  damped  modes  increase  with  drift 
velocity  in  the  range  15-125  m/s,  a feature  not  predicted  by  the 
linear  theory  of  the  gradient- drift  instability.  No  appreciable 
difference  is  seen  in  the  spectral  widths  of  the  upward  moving  waves 
as  compared  with  the  downward  moving  waves.  Moreover  the  width  to 
mean  ratios  of  the  computed  spectra  suggest  strongly  turbulent 
fluctuations  in  agreement  with  radar  observations  [Balsley  and  Farley, 
1971]. 

There  remains  some  noise-like  features  in  the  computed  spectral 
characteristics  of  the  short  wavelength  (X  < 10  m)  irregularities. 

This  is  to  be  expected  due  to  resolution  limitations.  In  the  models 
studied  a grid  mesh  of  64  x 64  points  is  used  which  allows  Fourier 
components  of  up  to  wavenumber  32  in  each  direction.  Components  with 
wavenumbers  greater  than  16,  therefore,  are  subject  to  possibly  serious 
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distortion  from  aliasing  errors.  The  highest  wavenumber  studied  in 
Model  1-3  was  14  (9.1  nt) . The  overall  trend  of  increasing  line- 
width  with  decreasing  wavelength  is  noted  in  the  decay  rates  of  the 
large  amplitude  short  wavelength  test  waves.  The  effects  on  the 
spectra  at  the  short  wavelengths  due  to  other  sources  of  error,  e.g., 
spatial  and  temporal  truncation  errors,  could  not  be  fully  investigated 
due  to  the  high  computational  cost  of  the  spectral  integration  times. 
However,  in  Model  3,  after  restarting  near  the  quasi-steady  state, 
several  short  wavelength  spectra  were  computed  after  halving  the 
time  step  thereby  reducing  the  time  truncation  error.  No  appreciable 
differences  were  found  in  spectral  features. 

Finally,  the  power  law  (1^  « k n,  n = 3-4)  of  the  turbulence  and 
the  scaling  of  the  frequency  broadening  with  wavenumber  k and  drift 
velocity  lend  support  to  the  strong  turbulence  theory  of  Sudan  and 
Keskinen  (1977)  which  predicts  k^  1^' ^ with  Ij.  k n,  n = 3.2. 
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APPENDIX  A: 

METHOD  OF  COMPUTING  TEMPORAL  POWER  SPECTRA 

Fourier  expanding  6n(x,t)  over  a s([uare  of  edge  L in  the  time 
interval  T 


6n(x,t)  = J’SnCk^w)  exp{i(k>x,  - cat)}  (A.l) 

k,(a 

6n(k,co)  = (2tt)  "^(L2T)  * f d2x  dt  Sn(x,t)  exp{-i(k*x  - cat)}  (A.2) 
*•  I '*»  w- 


where  x = (y,z),  jc^=  (k^,kz)  and  defining  the  normalized  temporal 
power  spectrum  Ij,  w 


< 


6n 

no 


d2k 


da)  I 


k,u) 


(A.  3) 


we  find  for  homogeneous  and  stationary  turbulence 


= lira 
L,T-*>° 


(A.  4) 


where  the  angle  brackets  denote  an  ensemble  average,  viz.,  the  average 
over  a large  number  of  identical  systems  of  which  the  system  under 
study  is  a member.  In  the  quasi-steady  state,  the  ensemble  average 
is  equivalent  to  a time  average  over  any  one  system  of  the  ensemble. 
This  direct  method  of  computing  the  temporal  power  spectra  is  more 
desirable  than  the  method  of  autocorrelations  due  to  computer  storage 
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limitation  and  the  relative  speed  of  the  fast  Fourier  transform  (FFT). 
Both  methods  can  he  shown  to  be  equivalent  [Blackman  and  Tukey,  1959], 
In  order  to  compute  I,  ^ in  practice  we  must  first  compute 
Snfk.uj)  which  in  discrete  form  can  be  written  for  fixed  k 

f+r  /*• 

N-l 

6n(lc,$.Af)  = ^ \ Sn(k,  jAt)exp{-2Trijfl,/N}  (A.  5) 

j=0 

where  N is  the  number  of  time  samples,  1 = 0,  1,  . . N - 1,  j = 0,  1, 
N - 1,  Af  = T,  T the  total  integration  time.  At  is  the  sampling 
period.  Power  spectral  estimates  from  successive  time  blocks  were 
averaged  to  reduce  the  variance.  The  sampling  frequency  (At)  * was 
chosen  to  be  always  greater  than  twice  the  frequencies  in  which 
there  was  measurable  power  in  order  to  minimize  chasing.  This  was 
accomplished  by  first  computing  the  highest  frequencies  in  the 
system.  In  our  case,  these  were  found  at  the  shortest  wavelengths 
(14  < k/kQ  > 10,  9 m < A < 12  m).  In  Madels  1-3,  the  maximum 
frequencies  corresponding  to  these  short  wavelengths  was  found,  by 
trial  and  error,  to  be  less  than  10  sec 
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SUMMARY  OF  M3DEL  PARAMETERS 


BEST  FIT  POWER  LAWS  AND  RELATIVE  COEFFICIENTS  FOR  M3  DEL  2 


Comparison  of  nonlinear  damping  as  derived  from  spectral  width,  test  wave  decay  rates. 
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F I CURE  CAPTIONS 

Figure  1.  Maximum  total  density  fluctuation  6n/n,  amplitude  of 

initially  groining  mode  q,  and  maxinum  horizontal  and  vertical 

turbulent  velocities  V and  V as  function  of  time  in  seconds. 

n v 

(a)  Model  1. 

(b)  Model  2. 

(c)  Model  3. 

Figure  2.  Contour  plots  of  density  fluctuation  <5n  for  Model  3 
(V^  = 125  msec  ^).  Contour  lines  represent  equal  increments  in 
<5n  between  maximum  (+  signs)  and  minimum  (-  signs).  Horizontal 
axis  is  y (direction  of  electron  drift)  and  vertical  axis  is 
z (vertical).  Scale  for  A^  ^ reads  one  half  of  6n/n  scale. 

(a)  t = 0,  6n/n  = .04,  Aj  Q = 6.14,  VH  max  = 100. 

(b)  t = 2,  <5n/n  = .06,  A1>0  = 9.22,  VH>max  = 173. 

Cc)  t « 4,  5n/n  - .14,  A1>0  ■ 13.61,  VH>MX  - 427. 

(d)  t - 12,  Sn/n  = .07,  Ax  0 - 5.81,  VH>max  - 251. 

(e)  t = 26,  6n/n  = .068,  A.  n = 4.83,  V„  i_„v  = 212. 

X)U  n,  nicix 

(f)  t = 34,  6n/n  = .06,  A1>()  = 4.53,  VH>mx  = 192. 

Figure  3.  Sample  time  averaged  Fourier  coefficients  I(k)  denoted  by 
closed  circles,  and  f 2 (^)  denoted  by  open  circles,  computed  near 
beginning  of  quasi-steady  states  of  Model  2. 

Figure  4.  Time  evolution  of  I(k)  vs.  k for  Model  2 (V^  = 100  msec  * 

(a)  t = 0 sec. 

(b)  t = 2 sec. 

(c)  t = 13  sec. 

(d)  t = 25  sec. 
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Figure  5.  Plot  of  I (k)  vs.  Vj  for  several  wavelengths  A in  the  quasi- 
steady  states  of  Models  1-3. 

Figure  6.  Comparison  between  computed  and  radar-obtained  power 
spectra  Ij.  ^ vs.  frequency  in  Hz.  Solid  curve  from  numerical  simu- 
lation; dotted  curve  from  radar  backscatter  observations  at  16.25 
NHz  (A  = 9.2  m)  [Balsley  and  Farley,  1971];  Aw/2it  is  the  difference 
of  the  frequencies  at  the  upper  and  lower  half-power  points. 

(a)  Radar  spectrum  received  from  60°  east  of  vertical;  simulation 
spectrum  taken  at  tan  1 (k^/k^)  = -30.2°. 

(b)  Radar  spectrum  received  from  45°  west  of  vertical;  simulation 

spectrum  taken  at  tan  1 (k /k  ) = 45°. 

y z 

(c)  Simulation  spectrum  taken  at  tan  1 (k  /k  ) =90°  and 

y z 

2 2 1/2 

2tt/ (ky,  + k^)  ' = 16  m.  Vertical  9 m spectra  have  not  been  studied 

with  radar  due  to  critical  frequency  limitations. 

Figure  7.  Computed  spectrum  frequencies  as  functions  of  k for 
Models  1-3.  The  spectrum  frequency  u)/2tt  (Hz)  is  the  average  of  the 
frequencies  at  the  upper  and  lower  half-power  points.  Solid-dashed 
line  is  frequency  from  linear  theory  at  = 125  msec  solid  line 
for  Vj  = 100  msec-1;  dashed  line  for  = 75  msec-1.  Arrow  marks 
critical  wavenumber  kc  from  linear  theory. 

(a)  0 = tan  1 (k  /k  ) = 90°;  linear  theory  predicts  nonpropagating 

y z 

vertical  modes. 

(b)  6 = 45°,  -45°. 

(c)  0 = 30°,  -30°. 

Figure  8.  Comparison  of  computed  phase  velocity,  w/k,  in  steady  state 
with  predicted  dependence  from  linear  theory  (w/k  ® cos  0)  for  Model 
2 (Vj  = 100  msec  1).  Solid  line  is  linear  theory  at  Vj  = 100  msec  . 


Figure  9.  Computed  spectrum  line  width  as  functions  of  k for  Models 
1-3.  Solid  line  is  linear  damping  rate.  Arrow  marks  critical  wave- 
number  kc  from  linear  theory. 

(a)  9 » 90°. 

(b)  0 = 45°,  -45°. 

(c)  0 = 30°,  -i0°. 

Figure  10.  Plot  of  Am/m  for  Models  1-3  at  0 = 45°. 

Figure  11.  Comparison  between  numerically  obtained  angle- averaged 
spectral  widths  and  predicted  angle- averaged  nonlinear  frequency 
broadening  vs.  k for  Itodels  1-3  Dots,  circles,  and  triangles 
show  Am  from  simulation  averaged  over  angles;  dotted  curve  gives 
Yj,  from  linear  theory;  solid  lines  show  (T£  + Yj.)  ; at  large  k, 

Yj,  is  comparable  to  Tj.  in  simulations;  numerical  factors  are  adjusted 
for  best  fit. 

Figure  12.  Plot  of  computed  angle  averaged  spectral  linewidth  Am 
vs.  Vj  at  several  different  wavelengths;  numerical  factors  are 
adjusted  for  best  fit. 

Figure  13.  Plot  of  Am  vs.  electrojet  drift  velocity  from  radar 
observations  at  45°  [Balsley  and  Farley,  1973]  for  X = 9 m. 

Figure  14.  Plot  of  decay  of  squared  amplitude  of  10  m test  wave  vs. 
time  in  seconds  launched  in  Model  2.  Decay  time  z is  taken  to  be 
time  at  which  anplitude  falls  to  e ^ of  initial  amplitude. 


400 


39 


t (sec) 


55 


3r- 


N 

X 

CVJ 


3 -I 


-2 


-3 


o ^d=75  m/sec 
• = 100  m/sec 

a Vd  = i25m/sec 


* 4 ° i * 


OA  o 

ft 

o 


o * Q 


i 


1 


1 


1 


1 


I 3 5 7 9 II 

128  42.7  25.6  18.3  14.2  11.6 


13  K/K0 
9.8  X(m) 


i 


w/2ir  (Hz)  w/2tt  (Hz) 


Acj/27t  (Hz) 


Agl>/2tt  (Hz)  Aw/2tt  (Hz) 


(zh)  M.z/mv 


< 


O 


o 


0 


9 A 


(27r)~  / 68  Acj  (Hz) 


1 


9 m 


V(m/s) 


t (sec) 


